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^ (57) Abstract: The invention concerns an adapted filter (5) produdng observations of a reference vector (r) of W successive samples 
^ of a signal formed by adapted filtering of a signal (Y) received through a radio propagation channel with a pilot spread code (Cr). 
^ A processing is performed in said adapted filter output to estimate the parameters representing the propagation channel, including 
^ specific vectors (VjJ of a self-correlating matrix of the reference vector and the specific associated values (A*). Said parameters can 

be used to perform coherent demodulation of the signals received through other logic channels sharing the same physical propagation 

channel. 



^ (57) Abr6ge: Un filtre adaptd (5) produit des observations d'un vecteur de rdfSrence (r) de W «;hantillon5 successifs d'un signal 
_ formd par filtrage adapts d*un signal (Y) re^u selon un canal de propagation radio avec un code d'etalemenl pilote (Cr). Un iraitement 

est effectue en sortie de ce filtre adapts pour estimer des paramfetres representant le canal de propagation, iiicluant des vecteurs 
Q propres (v) d'une riiatrice d*autocorr^ation du vecteur de r^f^rence et les valeurs propres associ^ (Xk). Ces param^tres peuvent 

eire utilises pour effectuer une demodulation coh^rente de signaux regus suivani d'auu-es canaux logiqucs partagearit le m6me canal 
^ physique de propagation. 
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DISPQSITIFS D' ESTIMATION DE CANAL RADIO ET DE DEMODULATION 



La presente invention concerne ies techniques d'estimation des canaux 
radio, ainsi qu'un demodulateur utilisant une telle technique, 

L'invention trouve application dans le domaine des 
5 radiocommunications a acces multiple a repartition par codes (CDMA, « Code- 
Division Multiple Access »), notamment dans Ies cas ou le temps d'un symbole 
d'information est sensiblement plus grand que la memoire du canal (c'est-a- 
dire qu'on peut negliger Tinterference inter-symboles). 

Dans un systeme CDMA, Ies symboles transmis, binaires (± 1 ) ou 
10 quaternaires (±1 ±j). sont multiplies par des codes d'etalement composes 
d'echantillons, appeles « chips », dont la cadence est superieure a celle des 
symboles. Des codes d'etalement orthogonaux ou quasi-orthogonaux sont 
alloues a differents canaux logiques partageant la meme frequence porteuse, 
afin de permettre a chaque recepteur de detecter la sequence de symboles qui 
15 lui est destinee. en multipliant le signal regu par le code d'etalement 
correspondant. 

Un avantage du systeme CDMA est qu'il permet d'exploiter des trajets 
multiples de propagation radio, en utilisant ce qu'on appelle un recepteur en 
rateau (« rake receiver »). 

20 Le recepteur en rateau effectue une demodulation coherente fondee 

sur une approximation de la reponse impulsionnelle du canal de propagation 
radio par une serie de pics, chaque pic apparaissant avec un retard 
correspondant au temps de propagation le long d'un trajet particulier et ayant 
une amplitude complexe correspondant a Tattenuation et au dephasage du 

25 signal le long de ce trajet (realisation instantanee du fading). En analysant 
plusieurs trajets de reception, c'est-a-dire en echantillonnant plusieurs fois la 
sortie d'un filtre adapte au code d'etalement du canal, avec des retards 
correspondant respectlvement a ces trajets, le recepteur en rateau obtient des 
estimations multiples des symboles transmis, qui sont combinees pour obtenir 

30 un gain en diversite. La combinaison peut notamment etre effectuee selon la 
methode dite MRC (« Maximum Ratio Combining »). qui pondere Ies 
differentes estimations en fonction des amplitudes complexes observees pour 
Ies differents trajets. 

Afin de permettre cette demodulation coherente, un canal pilote peut 

35 etre prevu pour I'estimation de la reponse impulsionnelle sous forme d'une 
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succession de pics. La reponse impulsionnelle est estimee au moyen d'un filtre 
adapte a un code d'etalement pilote avec lequel remetteur module une 
. sequence de symboles connus, par exemple des symboles a 1. Les positions 
des maxima de la sortie de ce filtre adapte donnent les retards utilises dans les 
5 doigts du recepteur en rateau, et les amplitudes complexes associees 
correspondent aux valeurs de ces maxima. Lorsque plusieurs trajets distincts 
ont des retards tres proches. par exemple avec des ecails inferieurs a la duree 
d'un chip, la recherche des maxima conduit generalement a ne prendre en 
compte qu'un seul de ces retards. 

10 Un but de la presente invention est d'exploiter une representation du 

canal de propagation radio qui soit moins approximative que dans un recepteur 
en rateau traditionnel, notamment dans le cas ou 11 existe plusieurs trajets de 
propagation ayant des retards proches. 

Selon I'invention, il est propose un dispositif d'estimation de 

15 parametres representant un canal de propagation radio, comprenant un filtre 
adapte pour produire des observations d'un vecteur de reference de W 
echantillons successifs d'un signal forme par filtrage adapte d'un signal regu 
selon le canal de propagation avec un code d'etalement pilote. et des moyens 
d'estimation de parametres representant le canal de propagation, incluant des 

20 vecteurs propres d'une matrice d'autocorrelation du vecteur de reference. Les 
parametres representant le canal de propagation peuvent en outre comprendre 
des valeurs propres de la matrice d'autocon-eiation du vecteur de reference, 
respectivement associees auxdits vecteurs propres avec lesquels elles sont 
estimees. 

25 Le dispositif emploie une description du canal de propagation sous 

forme de ses elements propres, ce qui procure une estimation plus riche que le 
modele de la ligne a retard (retards et amplitudes) pris en compte dans le 
recepteur en rateau traditionnel. Ainsi. chaque trajet de propagation peut etre 
pris en compte dans la reponse du canal au moyen d'une forme d'onde 
30 correspondent a un vecteur propre de la matrice d'autocorrelation et non 
simplement par un pic ponctuel associe a un retard. 

Dans un mode de realisation prefere, Testimation des vecteurs propres 
de la matrice d'autocorrelation comprend les etapes suivantes : 

- extraire du vecteur de reference M sous-vecteurs de reference composes 
35 chacun de S echantillons successifs dudit signal forme par filtrage 

adapte, M et S etant des nombres plus grands que 1, les M sous- 
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vecteurs de reference correspondant en moyenne a M plages 
energetiques du vecteur de reference ; 

- estimer les matrices d'autocorrelation respectives des M sous-vecteurs 
de reference ; . 

5 - detemniner des sous-vecteurs propres associes chacun a une valeur 
propre et a I'un des M sous-vecteurs de reference, cheque sous-vecteur 
propre etant, pour la valeur propre associee. un vecteur propre de la 
matrlce d'autocorrelation estimee du sous-vecteur de reference associe ; 

- former cheque vecteur propre estime de la matrlce d'autocorrelation du 
10 vecteur de reference en plagant un des sous-vecteurs propres aux 

positions des S echantillons du sous-vecteur de reference associe et des 
echantillons nuls aux W-S autres positions. 
Ceci pennet de realiser un compromis entre la finesse d'estimation du 
canal et la complexite des calculs requis, qui sont toutes deux des fonctions 
15 croissantes de S. ainsi que de M. Plus la capacite de calcul disponible est 
grande, plus les sous-vecteurs peuvent etre de taille importante, ce qui affine 
I'approximation de la reponse du canal. 

Le dispositif ci-dessus peut etre utilise cheque fbis qu'on a besoin de 
caracteriser un canal de propagation, en particulier pour effectuer une 
20 demodulation. 

Un autre aspect de la presente invention se rapporte ainsi a un 
demodulateur pour une reception selon un canal de propagation radio sur 
lequel sont definis plusieurs canaux logiques multiplexes incluant un canal 
pilote et au moins un canal de donnees, comprenant des moyens d'obtention 

25 de plusieurs observations d'un premier vecteur de W echantillons successifs 
d'un premier signal forme par filtrage adapte d'un signal re?u selon le canal de 
propagation au moyen d'un premier code d'etalement relatif au canal pilote. 
des moyens d'obtention d'un second vecteur de W echantillons successifs d'un 
second signal forme par filtrage adapte dudit signal re?u au moyen d'un second 

30 code d'etalement relatif au canal de donnees, des moyens d'estimation de 
parametres representant le canal de propagation, incluant des vecteurs 
propres d'une matrice d'autocon-elation du premier vecteur, et des moyens 
d'estimation de symboles pour estimer au moins un symbole d'information 
porte par le second signal sur la base des premier et second vecteurs et 

35 desdits parametres representant le canal de propagation. Ce demodulateur 
incorpore un dispositif d'estimation de canal tel qu'enonce precedemment. 
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Si le second signal porte des symboles d'infomiation binaires modules 
selon una modulation binaire a deplacement de phase, les moyens 
d'estimation de symboles peuvent etre agences pour estimer au moins un 
symbole d'infonnation binaire porte par le second signal par le signe de la 
5 partie reelle du nombre complexe : 




ou d est le nombre de vecteurs propres estimes de la matrice d'autocorrelation 
du premier vecteur, est un vecteur propre de la matrice d'autocorrelation du 

premier vecteur, associe a une valeur propre X^, Nq est une estimation de la 
10 puissance du bruit sur ie canal de propagation, r et x sont respectivement les 
premier et second vecteurs, p est un coefficient predetermine, {.)^ desrgne le 
transpose conjugue, et |.|| la norme d'un vecteur. 

Une vraisemblance dudit symbole d'information binaire estime peut 
etre determinee proportionnellement a la valeur absolue de la partie reelle du 
15 nombre complexe y. 

Si le second signal porte des symboles dlnfonmation quaternaires 
modules selon une modulation quaternaire a deplacement de phase, les 
moyens d'estimation de symboles peuvent etre agences pour estimer au moins 
un symbole d'information quaternaire porte par le second signal par deux bits 
20 respectivement definis par les signes des parties reelle et imaginaire du 
nombre complexe y defini par la relation (1). Des vraisemblances respectives 
de ces deux bits peuvent en outre etre determinees proportionnellement aux 
valeurs absolues des parties reelles et imaginalres du nombre complexe y. 

Le demodulateur peut egalement etre adapte pour traiter un signal regu 
25 selon le canal de propagation sous forme de plusieurs signaux captes au 
moyen d'antennes differentes, afin de procurer une diversite spatiale. 

D'autres particularites et avantages de la presente invention 
apparaitront dans la description ci-apres d'exemples de realisation non 
limitatifs, en reference aux dessins annexes, dans lesquels : 
30 - la figure 1 est un schema synoptique d'un exemple de recepteur de 
radiocommunication comportant un demodulateur selon Tinvention ; 
- la figure 2 est un schema d'un module d'estimation de canal du 
demodulateur de la figure 1 ; 
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- la figure 3 est un schema synoptique d'un module d'estimation de 
symboles du demodulateur de la figure 1 ; et 

- la figure 4 est un schema synoptique d'une unite de calcul du module de 
la figure 3. 

Le recepteur represents sur la figure 1 comporte un etage radio 1 qui 
regoit le signal radio capte par I'antenne 2 et le convertit en bande de base. Le 
signal en bande de base est numerise par un convertisseur analogique- 
numerique 3, puis fourni a un filtre de reception 4. Le filtre 4 assure un filtrage 
adapte a la mise en forme des signaux par Temetteur. II delivre un signal* 
numerique Y a raison d'un echantillon complexe par chip des codes 
d'etalement. 

En presence de trajets multiples de propagation (au nombre de d), la 
reponse h(t) du canal allant jusqu*au filtre de reception 4 peut s'ecrire : 

d-1 

h(t)= j;sk.a(t-tk) (2) 
k=0 

Les representent les differents retards des echos, les sont les realisations 
instantanees du fading variant au cours du temps, et a(t) designe Timpulsion 
formant la modulation. On suppose que le critere de Shannon est respecte 
lorsque les signaux sont echantillonnes a la frequence chip 1/t, et que la 
reponse du canal est nulle pour t > {W"1)T. Dans ce cas, la reponse du canal 
peut etre representee par un vecteur a W composantes complexes : 

d"-i 

= Z^k'^k (3) 
k=0 

Dans la relation (3), les y,^ sont les vecteurs propres de la matrice 
d'autocorrelation de h (soit R.y,^ = lij^.y,^ avec R = E(h.h^), ou E(.) designe 
I'esperance mathematique). et les aj^ sont des realisations instantanees de 
variables aleatoires gaussiennes complexes de variances egales aux valeurs 
propres [x^ de la matrice autocorrelation (E^|ak|^ j = |ik )• Cette decomposition 

de la matrice R donne des vecteurs propres orthogonaux (y|^*^.y|^. = 0 si k ;t k') 
et des valeurs propres reelles positives ou nulles puisque la -matrice est 
hermitienne. 



h{0) 
h(T) 

h((W^1)T) 
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Parmi les canaux logiques multiplexes sur ie canal radio, on a un canal 
de controle sur lequel sont transmis des symboles pilotes et un ou plusieurs 
canaux de donnees sur lequel sont transmis des symboles d'information. On 
suppose que la puissance d'emission sur Ie canal de controle est (3 fois 
5 superieure a la puissance d'emission d'un symbole d'information b sur un canal 
de donnees (cette derniere puissance sera consideree comme unitaire). Le 
coefficient p est connu a I'avance. Le signal est affecte d'un bruit blanc additif 
centre gaussien de variance Nq. 

Un filtre adapte 5 filtre le signal regu Y conformement au code 
10 d'etalement Cr du canal pilote. II en resulte pour chaque symbole pilote un 
vecteur r ayant W composantes complexes, donne par : 

d-1 

r = Vp2]°^kVK+n (4) 
k=0 



avec E^.n*^)= No-IJa matrice identite etant notee I . 

Un autre filtre adapte 6 filtre le signal re?u Y conformement au code 
15 d'etalement Cx du canal de donnees. 11 en resulte pour chaque symbole 
inconnu b un vecteur x ayant W composantes complexes, donne par : 

d-1 

x = b.2^ak.yk+0' (5) 



k=0 



avec E(]'.rf^)=No.l et E(i.n'^)= 0. 



On note X = 



le vecteur gaussien de taille 2W constitue par la 



20 concatenation des signaux regus sur les deux canaux. Sa densite de 
probabilite p(X) peut s'ecrire : 



1 



P(X) = .exp(^x'^K-'^ x) (6) 

7t2^.det(K) ^ - ^ 

ou K = , et det(.) designe le determinant. 

Le but du demodulateur est de maximiser par rapport a b la fonction de 
25 vraisemblance donnee par la probabilite conditionnelle p(b|x). 

Comme les realisations du bruit et du fading sont independantes. on a : 



) 
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d-i 



Nol = B 



2;(2^lkWk.w|;') 



ksO 



(7) 



en utilisant la decomposition en elements propresde la matrice avec 



B = 



#1 g 

0° b.l 



at Wk = 



1 



V2.|v^ 



. Par application du lemme d'inversion 



matricielle, on obtient : 



d-1 

!-Z 



k=0 No+2nk.||Bw„f 



(8) 



En prenant le logarithme de la probabilite (6), et en notant que le 
determinant det(K) est commun a tous les symboles possibles et que 

||BW|^| =(1+3)/2, le demodulateur doit minimiser la fonction suivante par 
rapport a b : 



d-1 

(i)(b) = y 



2^k 



ix'^BWk 



|y|;'.(#r+b.x|' 



(9) 



k=o No +(1+PK 

Lorsque les symboles sont des bits modules selon una modulation 
binaire a deplacement de phase (BPSK, « binary phase shift keying »), 
I'estimation.b d'un symbole b, qui minimise la fonction ^{h), est alors donnee 
par le signe de la partie reelle du nombre complexe : 



d-1 



(10) 



(b = sgn[Re(y)]), et la vraisemblance maximisee est proportionnelle a la valeur 
absolue de Re{y). Si un decodage a entrees souples (« soft input decoding ») 
est effectue en aval du demodulateur, la vraisemblance du bit estime b peut 
done etre mesuree par A = \Re{y]^ • 

Lorsque les symboles sont quaternaires et modules selon une 
modulation quaternaire a deplacement de phase (QPSK, « quaternary phase 
shift keying »), Testimation d'un symbole b = b-| + j.b2, qui minimise la fonction 

(|)(b), est formee par deux bits donnes par les signes respectifs des parties 
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reelle et imaginaire du nombre complexe y defini par la relation (10) 
(b = bi + j.b2 = sgn[Re(Y)]+j.sgn[lm(y)]), et les vraisemblances de ces deux bits 
sont respectivement proportionnelles aux valeurs absolues des parties reelle et 
imaginaire de y. Si un decodage a entrees souples est effectue en aval du 
demodulateur, les vraisemblances des deux bits estimes b-|,b2 peuvent done 
etre mesurees par A-| =|Re(y)[ et A2 =|lm(y)|. 

Lorsque le recepteur comporte deux antennes en entree de deux voies 
de traitement, chaque voie ayant deux filtres adaptes I'un au code Cr du canal 
pilote fournissant un vecteur r^^, r^, et Tautre au code Cx du canal de donnees 
fournissant un vecteur x^, x^, la fonction a minimiser (j)(b) peut tenir compte de 
la diversite spatiale, son expression n'etant plus donnee par (9) mais par : 



d-1 

(t)(b)=X 



k=O^No+(1 + |3)^j^ 



I No+(UP)^k 



avec les notations : r*= 



X = 



in 



K p*.K 



yk".(VPr"+bx")| 

J 

(11) 

, p etant le facteur de 



correlation entre les deux antennes (soit E(ak_n Ct,^£,)=p.|ik). les vecteurs y,^* 
et y,^- etant les vecteurs propres de E' respectivement associes aux valeurs 



propres et soit yj=-^ 



V2\p.yK/|p|J' 



1 



>/2\-pVk/|p|_ 



= 0 +|p|)^ik et ^k = (1 -|p|Vk pour 0 ^ k < d. 

Dans le cas de reponses impulsionnelles de petite dimension W, il est 
aise d'extraire les vecteurs propres de E par une methode algebrique classique 
travaillant directement sur une estimation de la matrice, par exemple la 
methode dite « Power Iteration » (voir Golub & Loane : « Matrix Computation »; 
editions Baltimore M.B.). 

Dans le cas de reponses impulsionnelles de plus grande dimension W, 
la complexite requise pour mettre a jour I'estimation de R rend cette methode 

irrealisable en pratique. On utilise alors une matrice R = Et.r") definie comme 
I'autocorrelation du vecteur r, liee a R par la relation : 
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R = p.R + NoJ (12) 
On verifie immediatement que R et R ont les memes vecteurs propres 
Yi^. Les valeurs propres de R sont de celles de R multipliees par |3 et 
decalees de Nq, soit = ^.^^ + Nq. Avec (10), on en deduit Texpression (1) du 

nombre complexe y servant a estimer les bits tf information. 

En general, la matrice R a une structure creuse et possede peu 
d'elements significatifs. Eile peut etre decomposee selon Tapproximation 



10 



15 



20 



suivante 



R = 



0 0 
0 R° 



R*^^ 0 



(13) 



c'est-a-dire sous forme de blocs de taille SxS. notes R' pour 0 < i < M, 
distribues le long de sa diagonale principale (1 < S « W). Les positions de ces 
M blocs le long de la diagonale de R correspondent aux zones les plus 
energetiques de la reponse impulsionnelle du canal. De meme, la matrice R 
peut etre approchee par : 



R = 



No- 
0 



Nn.l 



.M-1 



0 

Nq.I 



(14) 



en notant R = R' h-Nq-I I'observee bruitee de pour 0 < i < M. 



Cheque sous-matrice R admet une decomposition en vecteurs 

propres, identique a celle de R', qui peut etre obtenue rapidement par une 
methode classique puisque sa dimension S est petite devant celle de R. En 
notant d' la dimension du sous-espace signal relativement a la sous-matrice 
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M-1 



R , avec = d, cette decomposition de R s'ecrit : 
~ i=0 



j=0 ^ 



H 



(15) 



De meme que precedemment, les valeurs propres sont decalees de la 
puissance du bruit. Ainsi les vecteurs propres de R sont : 

0 



avec < 



fO<k<d-1 
0 < i < M, 0 < j < d* 



(16) 



et ses valeurs propres sont celles des M matrices R*, decalees de Nq : 
_(^k-No). (yNo) f0^k<d-l 



3 



avec 



0 < i < M. 0 < j < d' 



(17) 



Les sous-matrices R sont localisees la ou la reponse impulsionnelle 
concentre ses maxima d'energie. Pour avoir la position des M sous-matrices, 
10 supposees de taille fixe SxS, il suffit done d'avoir une estimation du profil en 
puissance du canal et d'en extraire les M portions de taille S les plus 
energetiques. Ce profil peut etre simplement calcule a partir de la moyenne du 
module carre de chaque echantillon du vecteur r Pour calculer les sous- 
matrices d'autocorrelation R*. on detemnine les M sous-vecteurs r' de 

15 dimension S correspdndant aux signaux regus sur le canal pilote aux positions 
precedemment determinees, et on calcule apres un certain nombre 

* 'H 

d'observations des moyennes des sous-matrices r'.r* . 

II se peut que certaines valeurs propres estimees ne soient pas 
significatives ou que Ton choisisse de ne retenir qu'un certain nombre fixe 
20 d'elements propres. On peut done, en fin de traitement, effectuer une selection 
des vecteurs propres associes aux plus puissantes valeurs propres. 

Le demodulateur de la figure 1 comprend un module d'estimation de 
canal 8, represents de fafon plus detaille sur la figure 2, qui calcule des 
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estimations des vecteurs propres Vj^ et des valeurs propres de la matrice 
d'autocorrelation R du vecteur r delivre par le filtre adapte 5, ainsi qu'une 
estimation de la puissance du bruit Nq. 

En reference a la figure 2, le module 8 comprend une unite 80 qui 
5 calcule des vecteurs de W composantes constituees par les modules carres 
des composantes des observations successives du vecteur r delivre par le filtre 
adapte 5. Une unite 81 calcule une moyenne courante du vecteur delivre par 
Tunite 80, par example sur une fenetre rectangulaire ou exponentielle. L'unite 
82 determine les M plages de S echantillons les plus energetiques du vecteur 
10 moyenne, afin de definir les positions des sous-vecteurs r' de taille S (0 < i < M) 
au sein des vecteurs r 

Une unite 83 extrait les M sous-vecteurs r* de chaque vecteur r, et en 
estime la matrice d'autocorrelation par un calcul de moyennes des produits 
terme a terme : 



Des unites de diagonalisation 84 extraient les valeurs propres des 

matrices R*. du moins celles de plus grands modules, et les vecteurs propres 
associes Vj*. Ce calcul peut etre effectue selon toute methode algebrique 

connue, par exemple « Power Iteration », qui ne pose pas de probleme de 
20 mise en ceuvre compte tenu de la reduction operee sur les dimensions des 

matrices a diagonaliser (SxS au lieu de WxW). L'unite 85 selectionne les d 

valeurs propres de plus grands modules parmi les valeurs propres Xj' 

fournies par les unites 84, et les vecteurs propres correspondants qu'elle 

construit a partir des Vj* conformement a la relation (16). 
25 L'unite 85 peut egalement determiner la puissance du bruit Nq, par 

exemple prise egale a la plus petite des valeurs propres Xj". D'autres methodes 

classiques d'estimation de Nq peuvent etre utilisees. 

Revenant a la figure 1. les parametres Vj^, (0 < i < d) et Nq estimes 

par le module 8 sont fournis au module 9 qui estime les symboles transmis a 
30 partir des vecteurs r et x delivres par les filtres adaptes 5, 6. 




(19) 



La figure 3 illustre une structure possible du module 9 dans le cas 
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d'une modulation QPSK. Un additionneur 91 calcule la quantite y selon la 
relation (1) a partir des d contributions a sommer, calculees par des unites 
respectives 90. 

Une telle unite 90, relative a un trajet k (0 ^ k < d), est illustree par la 
5 figure 4. Les produits scalaires (r".y,^) et (y,^^.x) sont respectivement calcules 
en 95 et 96. puis multiplies entre eux par le multiplieur complexe 97. Le bloc 98 

effectue le calcul du coefficient ^ — ^ — — (on suppose ici que les unites 

84 ont delivre les vecteurs propres sous forme orthonormee, soit 
Vjj =||Ynf =1), qu'un multiplieur complexe 99 multiplie par la sortie du 

10 multiplieur 97 pour adresser la contribution du trajet k a I'additionneur 91 . 

En sortie de I'additionneur 91, les unites 92,93 extraient 

respectivement les parties reelle et imaginaire de y, dont les signes b^, ^2 

estiment les bits constituant le symbole quatemaire, et dont les valeurs 

absolues foumissent eventuellement des mesures de vraisemblance A^, y\2. 
15 . Dans le cas d'une modulation BPSK, le module d'estimation de 

symboles 9 peut avoir la meme structure que sur la figure 3. sauf que Tunite 93 

n'est pas necessaire. 

Dans les dessins annexes, le demodulateur est represente sous forme 

synoptique afin de mettre en evidence les differentes etapes de caicuL En 
20 pratique, ce calcul peut etre effectue par un processeur de signal numerique 

programme en consequence, ainsi qu'il est usuel dans les demodulateurs 

numeriques. 

Les inventeurs ont teste les performances du demodulateur ci-dessus 
dans le cas du systeme CDMA retenu pour les reseaux cellulaires de troisieme 

25 generation (voir « UMTS : la generation des mobiles multimedia », P. Blanc et 
al., L'Echo des Recherches, No. 170, 4^ trimestre 1997 & trimestre 1998, 
pages 53-68), dans le cas de codes d'etalement de cadence 3,84 Mchip/s avec 
une cadence de symboles de 15 kbit/s (facteur d'etalement SF = 256). Dans le 
systeme considere, le canal de donnees est constitue par la partie reelle (voie 

30 I) d'un signal complexe module en QPSK, transmis par un terminal mobile vers 
une station de base. Le canal pilote correspondant fait partie d'un canal de 
controle constitue par la partie imaginaire (voie Q) du signal complexe module 
en QPSK. Le signal est traite par blocs de 666 p^s. a raison de 2560 chips par 
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bloc, ce qui represente 10 bits sur chaque voie avec SF = 256. a savoir 10 bits 
d'information sur la voie I, et 6 bits pilotes et 4 bits de commande sur la voie Q. 
Bien que la modulation radio soit QPSK, les symboles binaires transmis sur 
chaque voie 1, Q doivent etre consideres comme faisant t'objet d'une 
5 modulation effective de type BPSK. Dans ces simulations, la longueur de la 
reponse du canal etait d'environ 15 ^s (W = 56). Pour chaque bit estime sur la 
voie I a partir d'un signal issu d'une seule antenne, le vecteur r utilise dans la 
relation (1) etait celui resultant de la moyenne des 6 bits pilotes du bloc 
courant, soit p = 6. Le canal etait estime au moyen de d = 4 vecteurs propres 
'10 determines au moyen de M = 3 sous-matrices de dimension S = 6. Pour 
differentes configurations du canal, les simulations ont montre un gain de 1/4 a 
3/4 de decibel sur le rapport Eb/Ng (Eb etant I'energie par bit) pour un taux 

d'erreur binaire (BER) egal a 10'^ avant decodage correcteur d'erreurs. Au prix 
d'une augmentation moderee de la complexite du recepteur. ce gain est tres 
15 appreciable pour augmenter la portee radio du systeme. 
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REVENDICATIONS 

1. Dispositif d'estimation de parametres representant un canal de 
propagation radio, comprenant un filtre adapte (5) pour produire des 
observations d'un vecteur de reference (r) de W ectiantillons successifs d'un 
signal form6 par filtrage adapte d'un signal (Y) regu selon le canal de 
propagation avec un code d'etalement pilote (Cr), et des moyens (8) 
d'estimation de parametres representant le canal de propagation, incluant des 
vecteurs propres (y,^) d'une matrice d'autocon-elation du vecteur de reference. 

2. Dispositif selon la revendication 1, dans lequel les parametres 
representant le canal de propagation, estimes par les moyens d'estimation de 
parametres (8). incluent en outre des valeurs propres (X,^) de la matrice 

d'autocon-elation du vecteur de reference (r). respectivement associees auxdits 
vecteurs propres {v^^). 

3. Dispositif selon la revendication 1 ou 2, dans lequel les moyens 
15 d'estimation de parametres (8) comprennent : 

- des moyens (80-83) d'extraction de M sous-vecteurs de reference (r*) 
composes chacun de S echantillons successifs dudit signal forme par 
filtrage adapte, M et S etant des nombres plus grands que 1, les M sous- 
vecteurs de reference correspondant en moyenne a M plages 

20 energetiques du vecteur de reference (r) ; 

- des moyens (83) d'estimation des matrices d'autocon-elation respectives 

- i 

( R ) des M sous-vecteurs de reference ; 

- des moyens (84-85) de detemiination de sous-vecteurs propres (Vj') 
associes chacun a une valeur propre (Aj') et a I'un des M sous-vecteurs 
de reference (r"). chaque sous-vecteur propre etant. pour la valeur propre 
associee, un vecteur propre de la matrice d'autocorrelation estimee du 
sous-vecteur de reference associe ; 

- des moyens (85) de fomiation de chaque vecteur propre estime (y,^) de la 
matrice d'autocorrelation du vecteur de reference en plapant un des 
sous-vecteurs propres (vj') aux positions des S echantillons du sous- 
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vecteur de reference associe et des echantillons nuls aux W-S autres 
positions. 

4. Dispositif selon la revendication 3, dans lequel les moyens de 
determination des sous-vecteurs propres (85) choisissent lesdits sous-vecteurs 
propres (Vj'), parmi Tensemble des vecteurs propres des matrices 

^autocorrelation (R*) des M sous-vecteurs de reference tf), comme etant ceux 

qui sont respectivement associes aux valeurs propres de plus grands 
modules. 

5. Dispositif selon la revendication 4, dans lequel les moyens de 
determination des sous-vecteurs propres (85) determinent un nombre donne 
(d) de sous-vecteurs propres (vj'), associes aux d valeurs propres (;^j^) de plus 
grands modules. 

6. Dispositif selon I'une quelconque des revendications 1 a 5, dans 
lequel ledit signal re^u selon le canal de propagation est compose de plusieurs 
signaux captes au moyen d'antennes differentes. 

7. Demodulateur pour une reception selon un canal de propagation 
radio sur lequel sont definis plusieurs canaux logiques multiplexes incluant un 
canal pilote et au moins un canal de donnees, comprenant des moyens (5) 
d'obtention de plusieurs observations d'un premier vecteur (r) de W 
echantillons successifs d'un premier signal forme par filtrage adapte d'un signal 
(Y) re?u selon le canal de propagation au moyen d'un premier code 
d'etalement (Cr) relatif au canal pilote, des moyens (6) d'obtention d'un second 
vecteur (x) de W echantillons successifs d'un second signal forme par filtrage 
adapte dudit signal re?u au moyen d'un second code d'etalement (Cx) relatif au 
canal de donnees, des moyens (8) d'estimation de parametres representant le 
canal de propagation, incluant des vecteurs propres (vj^) d'une matrice 
d'autocorrelation du premier vecteur. et des moyens (9) d'estimation de 
symboles pour estimer au moins un symboie d'information porte par le second 
signal sur la base des premier et second vecteurs et desdits parametres 
representant le canal de propagation. 
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8. Demodulateur selon la revendication 7, dans lequel les moyens 
d'estimation de parametres (8) comprennent : 

- des moyens (80-83) d'extraction de M sous-vecteurs de reference (r') 
composes chacun de S echantillons successifs du premier signal. M et S 

5 etant des nombres plus grands que 1 , les M sous-vecteurs de reference 

correspondent en moyenne a M plages energetiques du premier vecteur 
(r); 

- des moyens (83) d'estimation des matrices d'autocorrelation respectives 
(r') des M sous-vecteurs de reference ; 

10 - des moyens (84-85) de determination de sous-vecteurs propres (Vj') 

associes chacun a une valeur propre {X'l) et a Tun des M sous-vecteurs 

de reference (r*), chaque sous-vecteur propre etant, pour la valeur propre 
associee, un vecteur propre de la matrice d'autocorrelation estimee du 
sous-vecteur de reference assocle ; 
15 - des moyens de formation de chaque vecteur propre estime {^) de la 

matrice d'autocorrelation du premier vecteur en plagant un des sous- 
vecteurs propres (Vj') aux positions des S echantillons du sous-vecteur de 
reference associe et des echantillons nuls aux W-S autres positions. 

9. Demodulateur selon la revendication 8. dans lequel les moyens de 
20 determination des sous-vecteurs propres (85) choisissent lesdits sous-vecteurs 

propres (Vj*), panni I'ensemble des vecteurs propres des matrices 

d'autocorrelation (r') des M sous-vecteurs de reference (H), comme etant ceux 
qui sont respectivement associes aux valeurs propres (Aj') de plus grands 
modules. 

25 10. Demodulateur selon la revendication 9, dans lequel les moyens (85) 
de determination des sous-vecteurs propres determinent un nombre donne (d) 
de sous-vecteurs propres (Vj'), associes aux d valeurs propres (Xj') de plus 
grands modules. 

11. Demodulateur selon Tune quelconque des revendications 7 a 10, 

30 dans lequel les parametres representant le canal de propagation, estimes par 
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les moyens d'estimation de parametres (8), incluent en outre des valeurs 
propres (^) de la matrice d'autocorreiation normalisee du premier vecteur, 
respectivement associees auxdits vecteurs propres (V|^). 

12. Demodulateur selon la revendication 11. dans lequel le second 
5 signal porte des symboles ^'information binaires (b) modules selon une 

modulation binaire a deplacement de phase, et dans lequel les moyens 
d'estimation de symboles (9) sont agences pour estimer au moins un symbole 
d'information binaire porte par le second signal par le signe (b) de la partie 

(?^k-No)/|kf /h V h \ 
reelle du nombre complexe > -7 — r — — ^fc Vk Ib^k xl, ou d est le nombre 

10 de vecteurs propres estimes de la matrice d'autocorrelation du premier vecteur, 
Vj^ est un vecteur propre de la matrice d'autocorrelation normalisee du premier 

vecteur, associe a une valeur propre X^, Nq est une estimation de la puissance 
du bruit sur le canal de propagation, r et x sont respectivement les premier et 
second vecteurs. p est un coefficient predetermine, {,p designe le transpose 
15 conjugue, et j|. | la norme d'un vecteur. 

13. Demodulateur selon la revendication 12, comprenant des moyens 
de determination d'une vraisemblance (A) dudit symbole d'information binaire 
estime (b), proportionnellement a la valeur absolue de la partie reelle dudit 
nombre complexe. 

20 14. Demodulateur selon la revendication 11, dans lequel le second 
signal porte des symboles d'information quatennaires (b) modules selon une 
modulation quaternaire a deplacement de phase, et dans lequel les moyens 
d'estimation de symboles (9) sont agences pour estimer au moins un symbole 
d'information quaternaire porte par le second signal par deux bits (b-), 82) 

25 respectivement definis par les signes des parties reelle et imaginaire du 

nombre complexe Y (^k-No)/|kk|| t^^Vkj^fx), oCj d est le nombre de 

vecteurs propres estimes de la matrice d'autocorrelation du premier vecteur, 
est un vecteur propre de la matrice d'autocorrelation normalisee du premier 
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vecteur. associe a une valeur propre X^, Ng est une estimation de la puissance 
du baiit sur le canal de propagation, r et x sont respectivement les premier et 
second vecteurs. p est un coefficient predetermine, {.)^ designe le transpose 
conjugue, et ||.|| la norma d'un vecteur. 

5 15. Demodulateur selon la revendication 14, comprenant des moyens 
de determination de vraisemblances respectives (A-|, A2) desdits bits 

(b-|. 62), proportionnellement aux valeurs absolues des parties reelle et 
imaginaire dudit nombre complexe. 

16. Demodulateur selon Tune quelconque des revendications 7 a 15. 
10 dans lequel ledit signal regu selon le canal de propagation est compose de 
plusieurs signaux captes au moyen d'antennes differentes. 
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CM clfcctue en sortie de ce filue adapts pour csiimer de$ param^res rcpr^sentant le canal dC pmpagaiion. inchiani des veaeurs 
' propies (v) d'unc manicc d'autocoirdaiion du vecteur de r^fdreiKC d les valeurs propies assodces (^). Ces param&tits peovem 

kue uiilisds pour eflectuer unc d6ncdulation coh^ente de signaux le^us suivani d^autrcs canaux logiques pana^cant le m£me canal 

physique de propagaiion. 



RADIO CHANNEL ESTIMATION AND DEMODULATION DEVICES 

The present invention relates to techniques for estimating radio 
channels, as well as to a demodulator using such a technique. 

The Invention finds application in the field of code-division multiple 
5 access (CDMA) radio communications, especially in cases where the time of 
an information symbol is substantially greater than the memory of the channel 
(i.e. inter-symbol inference can be neglected). 

In a CDMA system, the symbols sent, binary (±1) or quaternary 
(±1±j), are multiplied by spreading codes composed of samples, called 
Id "chips", whose rate is greater than that of the symbols. Orthogonal or quasi- 
orthogonal spreading codes are allotted to the various logical channels sharing 
the same carrier frequency, so as to allow each receiver to detect the sequence 
of symbols which is destined therefor, by multiplying the signal received by the 
corresponding spreading code. 

15 An advantage of the CDMA system is that it makes it possible to utilize 

radio propagation multipaths, by using what is referred to as a rake receiver. 

The rake receiver performs coherent demodulation based on an 
approximation of the impulse response of the radio propagation channel by a 
series of peaks, each peak appearing with a delay con-esponding to the 

20 propagation time along a particular path and having a complex amplitude 
conresponding to the attenuation and to the phase shift of the signal along this 
path (instantaneous realization of fading). By analyzing several reception paths, 
i.e. by sampling several times the output from a filter matched to the spreading 
code of the channel, with delays corresponding respectively to these paths, the 

25 rake receiver obtains multiple estimates of the symbols sent, which are 
combined to obtain a diversity gain. Combining can be performed in particular 
according to the so-called MRC (Maximum Ratio Combining) method, which 
weights the various estiniates as a function of the complex amplitude obsen/ed 
for the various paths. 

30 In order to allow this coherent demodulation, a pilot channel can be 
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provided for the estimation of the impulse response in the form of a succession 
of peaks. The impulse response is estimated by means of a filter matched to a 
pilot spreading code with which the transmitter modulates a sequence of known 
symbols, for example symbols set to 1 . The positions of the maxima of the 
5 output from this matched filter give the delays used in the fingers of the rake 
receiver, and the associated complex amplitudes correspond to the values of 
these maxima. When several distinct paths have very close delays, for example 
with discrepancies smaller than the duration of a chip, the search for the 
maxima generally leads to just one of these delays being taken into account. 

10 An object of the present invention is to utilize a representation of the 

radio propagation channel which is less approximate than in a traditional rake 
receiver, in particular in the case where several propagation paths having 
similar delays exist. 

According to the invention, there is proposed a device for estimating 
15 parameters representing a radio propagation channel, comprising a matched 
filter for producing observations of a reference vector of W successive samples 
of a signal formed by matched filtering a signal received along the propagation 
channel with a pilot spreading code, and means for estimating parameters 
representing the propagation channel, including eigenvectors of an 
20 autocorrelation matrix of the reference vector. The parameters representing the 
propagation channel may further Include eigenvalues of the autoconrelation 
matrix of the reference vector, which values are respectively associated with 
said eigenvectors with which they are estimated. 

The device employs a description of the propagation channel in the 
25 fonn of its eigenelements, thereby affording a richer estimate than the delay 
line model (delays and amplitudes) taken into account in the traditional rake 
. receiver. Thus, each propagation path may be taken Into account in the 
response of the channel by means of a waveform corresponding to an 
eigenvector of the autocorrelation matrix and not simply through a pointwise 
30 peak associated with a delay. 

In a preferred embodiment, the estimation of the eigenvectors of the 
autocorrelation matrix comprises the following steps: 
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- extracting from the reference vector M reference subvectors each 
composed of S successive samples of said signal formed by matched 
filtering, M and S being numbers greater than 1 , the M reference 
subvectors corresponding on average to M energy ranges of the 

5 reference vector; 

- estimating the respective autocorrelation matrices of the M reference 
subvectors; 

- determining the eigensubvectors each associated with an eigenvalue and 
with one of the M reference subvectors, each eigensubvector being, for 

10 the associated eigenvalue, an eigenvector of the estimated 

autocorrelation matrix of the associated reference subvector; 

- forming each estimated eigenvector of the autocorrelation matrix of the 
reference vector by placing one of the eigensubvectors at the positions of 
the S samples of the associated reference subvector and zero samples at 

1 5 the W-S other positions. 

This allows a compromise to be achieved between the fineness of 
estimation of the channel and the complexity of the calculations required, both 
of which are increasing functions of S, and of M. The larger the computational 
capacity available, the bigger may be the size of the subvectors, thereby 
20 sharpening the approximation of the response of the channel. 

The above device can be used whenever it is necessary to 
characterize a propagation channel, in particular to perform demodulation. 

Another aspect of the present invention thus pertains to a demodulator 
for reception along a radio propagation channel on which are defined a plurality 

25 of multiplexed logical channels including a pilot channel and at least one data 
channel, comprising means for obtaining a plurality of observations of a first 
vector of W successive samples of a first signal fomned by matched filtering a 
signal received along the propagation channel by means of a first spreading 
code relating to the pilot channel, means for obtaining a second vector of W 

30 successive samples of a second signal formed by matched filtering said signal 
received by means of a second spreading code relating to the data channel, 
means for estimating parameters representing the propagation channel, 
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including eigenvectors of an autoconrelation matrix of the first vector, and 
symbol estimation means to estimate at least one infomriation symbol earned 
by the second signal on the basis of the first and second vectors and of said 
parameters representing the propagation channel. This demodulator 
5 incorporates a channel estimation device as stated hereinabove. 

If the second signal canies binary infomiation symbols modulated by 
binary phase shift keying, the symbol estimation means may be arranged to 
estimate at least one binary infonnation symbol carried by the second signal by 
the sign of the real part of the complex number: 



10 




(1) 



where d is the number of estimated eigenvectors of the autocorrelation matrix 
of the first vector, is an eigenvector of the normalized autocorrelation matrix 

of the first vector, associated with an eigenvalue X^^, Ng is an estimate of the 
noise power on the propagation channel, r and x are respectively the first and 
15 second vectors, p is a predetermined coefficient, (.)^ designates the conjugate 
transpose, and || . || the nomi of a vector. 

A likelihood of said estimated binary information symbol can be 
detemiined proportionally to the absolute value of the real part of the complex 

ft 

number y. 

20 If the second signal canies quaternary infonnation symbols modulated 

by quartemary phase keying, the symbol estimation means may be arranged to 
estimate at least one quaternary infomiatipn symbol carried by the second 
signal by two bits respectively defined by the signs of the real and imaginary 
parts of the complex number y defined by relation (1). Respective likelihoods of 

25 these two bits can further be detemiined proportionally to the absolute values 
of the real and imaginary parts of the complex number y. 

The demodulator can also be adapted to process a signal received 
along the propagation channel in the fomn of a plurality of signals picked up by 
means of different antennas, so as to afford space diversity. 
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Other features and advantages of the present invention will become 
apparent in the following description of non-limiting exemplary embodiments, 
with reference to the appended drawings, in which: 

- figure 1 is a block diagram of an exemplary radio communication receiver 
5 comprising a demodulator according to the invention^ 

- figure 2 is a diagram of a channel estimation module of the demodulator 
of figure 1; 

- figure 3 is a block diagram of a symbol estimation module of the 
demodulator of figure 1; and 

.10 - figure 4 is a block diagram of a calculation unit of the module of figure 3. 

The receiver shown in figure 1 comprises a radio stage 1 which 
receives the radio signal picked up by the antenna 2 and converts it into 
baseband. The baseband signal is digitized by an analogue/digital converter 3, 
then supplied to a reception filter 4. The filter 4 carries out filtering matched to 
15 the shaping of the signals by the transmitter. It delivers a digital signal Y as one 
complex sample per chip of the spreading codes. 

In the presence of multiple propagation paths (d in number), the 
response h(t) of the channel up to the reception filter 4 can be written: 

d-1 

h(t)=5]sk.a(t-xk) (2) 
k=0 

20 The Tk represent the various delays of the echoes, the s^ are the instantaneous 

realizations of the time varying fading and a(t) designates the shaping pulse of 
the modulation. It is assumed that the Shannon criterion is complied with when 
the signals are sampled at the chip frequency 1/T, and that the response of the 
channel is zero for t > (W-1)T. In this case, the response of the channel can be 
25 represented by a vector with W complex components: 

d-1 

k=0 . 

In relation (3), the y^ are the eigenvectors of the autocorrelation matrix 



h = 



h(0) 
h(T) 

h((W-1)T) 
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of h (i.e. E-Vk = Hk-v^ with R = E(h.h"J. where E(.) designates the mathematical 
expectation) and the are instantaneous realizations of complex Gaussian 
random variables whose variances are equal to the eigenvalues of the 

autocorrelation matrix (Eflakf ] = ^k)- This decomposition of the matrix B 



5 gives orthogonal eigenvectors (Vk^.Vk.. = 0 if k 9t k') and real positive or zero 
eigenvalues since the matrix is hermitian. 

Among the logical channels multiplexed on the radio channel, there is a 
control channel on which pilot symbols are sent and one or more data channels 
on which information symbols are sent. It is assumed that the transmission 
10 power on the monitoring channel is p times greater than the transmission 
power for an information symbol b on a data channel (the latter will be regarded 
as unitary). The coefficient p is known in advance. The signal is impaired by 
centered additive Gaussian white noise of variance Nq. 



15 spreading code Cr of the pilot channel. This results, for each pilot symbol, in a 
vector r having W complex components, and given by: 




A matched filter 5 filters the signal received Y in accordance with the 




(4) 



Another matched filter 6 filters the signal received Y in accordance with 
20 the spreading code Cx of the data channel. This results, for each unknown 
symbol b. In a vector x having W complex components and given by: 

d-l 




k=0 




7- 



10 



15 



We note by X = 



the Gaussian vector of size 2W constructed by 



concatenating the signals received on the two channels. Its probability density 
p(X) can be written: 

1 



P(X) = 



7r2W.det(K) 



,exp(-X"K-''x) 



(6) 



where K = , and det(.) designates the detenminant. 

The aim of the demodulator is to maximize with respect to b the 
likelihood function given by the conditional probability p(b | X). 



Since the realizations of the noise and of the fading are independent, 



we have: 



K = 



p.R b*VpR 
bVp.R R 



+ Nol=B, 



d-1/ X 



k=0 



B^+NqI 



(7) 



using the decomposition of the matrix B into eigenelements, with 



B = 



#1 0 
0^ b.l 



and Wk = 



1 



. By applying the matrix inversion 



lemma, we obtain: 



Nn 



d-1 



k=0 No+2nk-||BWk|| 



(8) 



Taking the logarithm of the probability (6), and noting that the 
determinant det(K) is common to all the possible symbols and that 

l^wj^l =(l+p)/2, the demodulator must minimize the following function with 

respect to b: 
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^(b) =y — — x"bw/ 

= y.Ji^M_|vH(vp,^b.xf (9) 

When the symbols are each modulated according to binary phase shift 
keying (BPSK), the estimate b of a symbol b, which minimizes the function 
5 <|»(b), is then given by the sign of the real part of the complex number: 

|H Yh,) 

(b = sgn[Re(y)]), and the maximized likelihood is proportional to the absolute 
value of Re(Y). If soft input decoding is performed downstream of the 
demodulator, the likelihood of the estimated bit b can therefore be measured 
10 byA = |Re(y)|. 

When the symbols are quaternary and modulated according to 
quatemary phase shift keying (QPSK), the estimate of a symbol b = b^ + j.b2, 

which minimizes the function <|)(b), is formed by two bits given by the respective 
signs of the real and imaginary parts of the complex number y defined by 
15 relation (10) (b = bi +j.b2 =sgn[Re(y)]+j.sgn[lm(y)]). and the likelihoods of 
these two bits are respectively proportional to the absolute values of the real 
and imaginary parts of y. If a soft Input decoding is performed downstream of 
the demodulator, the likelihoods of the two estimated bits 61,62 may therefore 
be measured by A-i = |Re(y)| et A2 = |lni(y)| . 

20 When the receiver comprises two antennas at the input of two 

processing pathways, each pathway having two matched filters, one to the 
code Cr of the pilot channel supplying a vector r^, r^, and the other to the code 

Cx of the data channel supplying a vector x^, x^j, the function to be minimized 

<t>(b) can take account of space diversity, its expression no longer being given 
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by (9) but by: 



*(b)=x; 



k=c(No+(1 + PK 



yJ".(VPf+bx') 



/\hf 



with the notation: r'= 







fx ^ 




"K 


p*K' 




. x'= 


' 5- 






a., 








pX 


K 



No+(1 + PK' " , 

(11) 

, p being the correlation 



factor for the two antennas (i.e. E(ak n «k (j)=P4^k)« vectors V|^* and 
being the eigenvectors of B' respectively associated with the eigenvalues 

• ^^k=(l + |p|K and 



±( 1 




V2lpyk/|p|, 


■ - Vfi-P-VK/IPlJ 



and i.e. yj^ = 

Hk =(l-|p|)^kforO<k<d. 

In the case of impulse responses of small dimension W, It is easy to 
10 extract the eigenvectors of R by a conventional algebraic method working 

directly on an estimate of the matrix, for example the so-called Tower Iteration" 
method (see Golub & Loane: "Matrix Computation", published by Baltimore 
M.B.). 

In the case of impulse responses of larger dimension W, the complexity 
15 required to update the estimate of 6 makes this method unrealizable in 

practice. A matrix R = E^.r'^) is then used, defined as the autocorrelation of the 

vector r, and related to B through the relation; 

R = j3.R + No.l (12) 

It is readily apparent that R and B have the same eigenvectors Vk- The 
20 eigenvalues Xk of R are those of B multiplied by p and shifted by Nq, i.e. 
A^k = P-Mk Nq. With (10), we deduce therefrom the expression (1) for the 
complex number y serving to estimate the inforrhation bits. 

In general, the matrix B has a sparse structure and possesses few 
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significant elements, it can be decomposed according to the following 
approximation: 



R = 



0 0 
0 R° 



0 • 



r"^-^ 0 

0 0 



(13) 



i.e. in the fomi of blocks of size S^S, denoted B' for 0 < i < M, distributed along 
its main diagonal (1 < S « W). The positions of these M blocks along the 
diagonal of R corespond to the zones of greatest energy of the impulse 
response of the channel. Likewise, the matrix R can be approximated by: 



R = 



Nn.l 0 



0 



No. I 



;M-1 



0 Nq.I 



(14) 



on denoting r' =R' +No.I the noisy observation of R' for 0 < i < M. 



10 
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Each submatrix R admits an eigenvector decomposition identical to 

that of R', which can be obtained rapidly by a conventional method since its 
dimension S is small compared with that of R. Denoting by d' the dimension of 



M-1 



the signal subspace relative to the submatrix; R, with ^d' = d, this 

i=o 



decomposition of R can be written: 



j=0 ^ 



(15) 
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Just as previously, the eigenvalues are shifted by the noise power. 
Thus the eigenvectors of B are: 

0 



with 



0<k<d-1 



0<i<M, 0<j<d' 



(16) 



and its eigenvalues are those of the M matrices R , shifted by Nq : 



l0<i<M, 0<j<d' 



(17) 



The submatrices R are located precisely where the impulse response 

concentrates its energy maxima. To get the position of the M submatrices, 
which are assumed to be of fixed size SxS, it is therefore sufficient to get an 
estimate of the power profile of the channel and to extract therefrom the M 
10 portions of size S which have the greatest energy. This profile can be simply 
calculated from the average of the squared modulus of each sample of the 

vector r. To calculate the autoconrelation submatrices r' , we detenmine the M 

subvectors H of dimension S corresponding to the signals received on the pilot 
channel at the previously detennined positions, and we calculate after a certain 

■ 'H 

15 numberof observations averages of the submatrices r*.r' . 

It may happen that certain estimated eigenvalues are not significant or 
that we chose to retain only a certain fixed number of eigenelements. We can 
therefore, at the end of processing, make a selection of eigenvectors 
associated with the most powerful eigenvalues. 

20 The demodulator of figure 1 comprises a channel estimation module 8, 

represented in greater detail in figure 2, which calculates estimates of the 
eigenvectors V|^ and of the eigenvalues of the autocorrelation matrix R of the 
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vector r delivered by the matched filter 5, as well as an estimate of the noise 
power Nq. 

With reference to figure 2, the module 8 comprises a unit 80 which 
calculates vectors of W components consisting of the squared moduli of the 

5 components of the successive observations of the vector r delivered by the 
matched filter 5. A unit 81 calculates a running average of the vector delivered 
by the unit 80, for example over a rectangular or exponential window. The unit 
82 detemnines the M ranges of S samples of greatest energy of the average 
vector, so as to define the positions of the subvectors H of size S (0 < i < M) 

10 within the vectors r. 

A unit 83 extracts the M subvectors r* of each vector r, and estimates 
the autocorrelation matrix thereof by calculating averages of the tenn-by-term 
products: 



r'=e 



f .H^ 



r'.r' 



(19) 



15 Diagonalization units 84 extract the eigenvalues X,j' of the matrices R', 

at least those of greatest moduli, and the associated eigenvectors y|'. This 

calculation can be performed according to any known algebraic method, for 
example "Power Iteration", which poses no implementational problem in view of 
the reduction effected on the dimensions of the matrices to be diagonalized 
20 (SxS instead of WxW). Jhe unit 85 selects the d eigenvalues of largest moduli 
from among the eigenvalues supplied by the unit 84, and the 

corresponding eigenvectors which it constructs from the Vj' in accordance 
with relation (16). 

The unit 85 can also determine the noise power Nq, for example taken 
25 equal to the smallest of the eigenvalues Xj'. Other conventional methods of 
estimating Nq can be used. 

Retuming to figure 1, the parameters y|^,A,|^ (0 < i< d) and Nq estimated 
by the module 8 are supplied to the module 9 which estimates the symbols sent 
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on the basis of the vectors r and x delivered by the matched filters 5, 6. 

Figure 3 illustrates a possible structure of the module 9 in the case of 
QPSK modulation. An adder 91 calculates the quantity y according to relation 
(1) froni the d contributions to be summed, calculated by respective units 90. 

Such a unit 90, relating to a path k (0 < k < d), is illustrated by figure 4. 
The scalar products (H^.v^) and (Vk'^.x) are calculated at 95 and 96 
respectively, then multiplied together by the complex multiplier 97. The block 98 

perfomis the calculation of the coefficient , . is assumed here 

that the units 84 have delivered the eigenvectors in orthonormal forni, i.e. 

2 o 

V • = ||vk|| = 1 ), which a complex multiplier 99 multiplies by the output of the 
multiplier 97 so as to address the contribution of path k to the adder 91. 

At the output of the adder 91, the units 92, 93 extract the real and 
imaginary parts respectively of y, whose signs bi, estimate the bits 
constituting the quaternary symbol, and whose absolute values may supply 
15 likelihood measurements , A2. 

In the case of BPSK modulation, the symbols estimation module 9 can 
have the same structure as in figure 3, except that the unit 93 is not necessary. 

In the appended drawings, the demodulator is represented in 
schematic fomi so as to highlight the various computational steps, in practice, 
20 this computation can be performed by a digital signal processor programmed 
accordingly, as is standard practice in digital demodulators. 

The inventors have tested the performance of the above demodulator 
in the case of the CDMA system adopted for third-generation cellular networks 
(see "UMTS: la generation des mobiles multimedia" [the generation of the 
25 multimedia mobiles], P. Blanc et al., L'Echo des Recherches, No. 170, 4**^ 
quarter 1997 & 1^* quarter 1998, pages 53-68), In the case of spreading codes 
of rate 3.84 Mchip/s with a symbol rate of 15 kbit/s (spreading factor SF = 256). 
In the system considered, the data channel consists of the real part (I pathway) 



■) 
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of a QPSK-modulated complex signal sent by a mobile terminal to a base 
station. The con-esponding pilot channel fomns part of a control channel 
consisting of the imaginary part (Q pathway) of the QPSK-modulated complex 
signal. The signal is processed in blocl<s of 666 ps, at a rate of 2560 chips per 

5 block, this representing 10 bits on each pathway with SF = 256, namely 10 
infomiatlon bits on the I pathway and 6 pilot bits and 4 control bits on the Q 
pathway. Although the radio modulation is QPSK, the binary symbols sent on 
each pathway I, Q should be regarded as forming the subject of an effective 
modulation of BPSK type. In these simulations, the length of the response of 

10 the channel was about 15 ps (W = 56). For each bit estimated on the I pathway 
from a signal emanating from a single antenna, the vector r used in relation (1) 
was that resulting from the average of the 6 pilot bits of the cunrent block, i.e. p 
= 6. The channel was estimated by means of d = 4 eigenvectors determined by 
means of M = 3 submatrices of dimension S = 6. For various configurations of 

15 the channel, the simulations have shown a gain of 1/4 to 3/4 decibels in the 
Eb/No ratio (Eb being the energy per bit) for a binary enror rate (BER) equal to 

10"2 before eror con-ecting decoding. At the cost of a moderate increase in the 
complexity of the receiver, this gain is very appreciable In increasing the radio 
range of the system. 
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CLAIMS 



1 . A device for estimating parameters representing a radio propagation 
channel, comprising a matched filter (5) for producing observations of a 
reference vector (r) of W successive samples of a signal formed by matched 

5 filtering a signal (Y) received along the propagation channel with a pilot 
spreading code (Cr), and means (8) for estimating parameters representing the 
propagation channel, including eigenvectors (y|^) of an autocorrelation matrix of 

the reference vector. 

2. The device as claimed in claim 1 , wherein the parameters 
10 representing the propagation channel, estimated by the parameter estimation 

means (8), further include eigenvalues {7^) of the autocorrelation matrix of the 
reference vector (r), respectively associated with said eigenvectors (^). 

3. The device as claimed in claim 1 or 2, wherein the parameter 
estimation means (8) comprise: 

15 - means (80-83) for extracting M reference subvectors (r') each composed 
of S successive samples of said signal formed by matched filtering, M 
and S being numbers greater than 1, whereby the M reference 
subvectors correspond on average to M energy ranges of the reference 
vector (r); 

20 - means (83) for estimating the respective autoconrelation matrices (r') of 
the M reference subvectors; 

- means (84-85) for determining eigensubvectors (Vj') each associated with 

an eigenvalue (Xj') and with one of the M reference subvectors (r"), 

whereby each eigensubvector is, for the associated eigenvalue, an 
25 eigenvector of the estimated autocorrelation matrix of the associated 

reference subvector; 

- means (85) for fomiing each estimated eigenvector (Vj^) of the 

autocon-elation matrix of the reference vector by placing one of the 
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eigensubvectors (Vj') at the positions of the S samples of the associated 
reference subvector and zero samples at the W-S other positions. 

4. The device as claimed in claim 3, wherein the means for determining 
the eigensubvectors (85) choose said eigensubvectors (Vj') from the set of 

^ i 

5 eigenvectors of the autocorrelation matrices (R ) of the M reference subvectors 
(H), as being those which are respectively associated with the eigenvalues (Xj') 
of largest moduli. 

5. The device as claimed in claim 4, wherein the means for detemiining 
the eigensubvectors (85) determine a given number (d) of eigensubvectors (yj'), 

10 associated with the d eigenvalues (Xj) of largest moduli. 

6. The device as claimed in any one of claims 1 to 5, wherein said 
signal received along the propagation channel is composed of a plurality of 
signals picked up by means of different antennas. 

7. A demodulator for reception along a radio propagation channel on 
15 which are defined a plurality of multiplexed logical channels including a pilot 

channel and at least one data channel, comprising means (5) for obtaining a 
plurality of observations of a first vector (r) of W successive samples of a first 
signal formed by matched filtering a signal (Y) received along the propagation 
channel by means of a first spreading code (Cr) relating to the pilot channel, 
20 means (6) for obtaining a second vector (x) of W successive samples of a 
second signal fomried by matched filtering said signal received by means of a 
second spreading code (Cx) relating to the data channel, means (8) for 
estimating parameters representing the propagation channel, including 
eigenvectors (y,^) of an autocorrelation matrix of the first vector, and symbol 

25 estimation means (9) to estimate at least one infomiation symbol carried by the 
second signal on the basis of the first and second vectors and of said 
parameters representing the propagation channel. 
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8. The demodulator as claimed in claim 7, wherein the parameter 

estimation means (8) comprise: 

- means (80-83) for extracting M reference subvectors (r*) each composed 
of S successive samples of the first signal. M and S being numbers 

5 greater than 1, whereby the M reference subvectors con-espond on 

average to M energy ranges of the first vector (r); 

- j 

- means (83) for estimating the respective autocorrelation matrices (R ) of 
the M reference subvectors; 

- means (84-85) for determining eigensubvectors (vj') each associated with 

10 an eigenvalue (Xj) and with one of the M reference subvectors (r*), 

whereby each eigensubvector is, for the associated eigenvalue, an 
eigenvector of the estimated autocorrelation matrix of the associated 
reference subvector; 

- means for fomriing each estimated eigenvector (y,^) of the autoconrelation 



15 matrix of the first vector by placing one of the eigensubvectors (Vj') at the 

positions of the S samples of the associated reference subvector and 
^ zero samples at the W-S other positions. 

9. The demodulator as claimed in claim 8, wherein the means for 



determining the eigensubvectors (85) choose said eigensubvectors (Vj') from 

20 the set of eigenvectors of the autocorrelation matrices (r' ) of the M reference 
subvectors (H), as being those which are respectively associated with the 
eigenvalues (Xj') of largest moduli. 

10. The demodulator as claimed in claim 9, wherein the means (85) for 
detemriining the eigensubvectors detennine a given number (d) of 

25 eigensubvectors (\^'), associated with the d eigenvalues {X *) of largest moduli. 

11. The demodulator as claimed in any one of claims 7 to 10, wherein 
the parameters representing the propagation channel, estimated by the 
parameter estimation means (8), further include eigenvalues {\) of the 



•) 

f 
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normallzed autocorrelation matrix of the first vector, respectively associated 
with said eigenvectors (y,^). 

12. The demodulator as claimed in claim 11, wherein the second signal 
carries binary information symbols (b) modulated by binary phase shift keying, 

5 and wherein the symbol estimation means (9) are arranged to estimate at least 
one binary information symbol canied by the second signal by the sign (b) of 

the real part of the complex number ^ where d is 

the number of estimated eigenvectors of the autoconrelation matrix of the first 
vector, Vi^ is an eigenvector of the nomrialized autocon-elation matrix of the first 

10 vector, associated with an eigenvalue X^^, Nq is an estimate of the noise power 
on the propagation channel, r and x are respectively the first and second 
vectors, p is a predetemnined coefficient, designates the conjugate 
transpose, and || . || the nonn of a vector. 

13. The demodulator as claimed in claim 12, cohiprising means for 
15 detemnining a likelihood {k) of said estimated binary information symbol (b), 

proportionally to the absolute value of the real part of said complex number. 

14. The demodulator as claimed in claim 11, wherein the second signal 
canies quatemary information symbols (b) modulated by quarternary phase 
keying, and wherein the symbol estimation means (9) are arranged to estimate 

20 at least one quatemary information symbol earned by the second signal by two 
bits (61,62) respectively defined by the signs of the real and imaginary parts of 

the complex number ^^^^^^i^^^-^^^^^' ^ '® number of 

k — 0 

estimated eigenvectors of the autocorrelation matrix of the first vector, y,^ is an 
eigenvector of the normalized autoconrelation matrix of the first vector, 
25 associated with an eigenvalue X^^, Nq is an estimate of the noise power on the 
propagation channel, r and x are respectively the first and second vectors, p is 
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a predetermined coefficient, designates the conjugate transpose, and || . || 
the norm of a vector. 

15. The demodulator as claimed in claim 14. comprising means for 
determining respective likelihoods (X^, Xg) of said bits (b-|, 62), proportionally 

5 to the absolute values of the real and imaginary parts of said complex number. 

16. The demodulator as claimed in any one of claims 7 to 15, wherein 
said signal received along the propagation channel is composed of a plurality of 
signals picked up by means of different antennas. 
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